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Abstract. The power spectral density in Fourier frequency domain, and the 
different variants of the Allan deviation (ADEV) in dependence on the averaging 
time are well established tools to analyse the fluctuation properties and frequency 
instability of an oscillatory signal. It is often supposed that the statistical 
uncertainty of a measured average frequency is given by the ADEV at a well 
considered averaging time. However, this approach requires further mathematical 
justification and refinement, which has already been done regarding the original 
ADEV for certain noise types. Here we provide the necessary background to 
use the modified Allan deviation (modADEV) and other two-sample deviations 
to determine the uncertainty of weighted frequency averages. The type of two- 
sample deviation used to determine the uncertainty depends on the method 
used for determination of the average. We And that the modADEV, which is 
connected with A-weighted averaging, and the two sample deviation associated 
to a linear phase regression weighting (parADEV) are in particular advantageous 
for measurements, in which white phase noise is dominating. Furthermore, we 
derive a procedure how to minimize the uncertainty of a measurement for a 
typical combination of white phase and frequency noise by adaptive averaging 
of the data set with different weighting functions. Finally, some aspects of the 
theoretical considerations for real-world frequency measurement equipment are 
discussed. 
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1. Introduction 

The central outcome of a frequency measurement is 
an average frequency and its associated uncertainty. 
The determination of the uncertainty depends on 
the method to determine the average frequency. 
The topic of this paper is the stringent derivation 
of the relations between averaging methods and 
associated uncertainty, and the identification of 
averaging methods, which are advantageous in the 
sense of approaching the tightest uncertainty for given 
noise types. 

The average frequency value determined in a 
single measurement can only be an estimate of the 
true frequency, because even if there are no unknown 
systematic offsets, we only have limited knowledge 
about the actual value due to noise which always 
contributes to a measurement. Often, the term 
“statistical uncertainty” is used with this meaning, and 
in the following, we will use the term “uncertainty” or 
“tt” as a short form. 

Traditionally, frequency has been most often 
measured with dead-time free so called Il-counters. 
The most natural measure of frequency instability 
determined from Il-counted frequency values is the 
Allan deviation (ADEV) [1, 2]. Besides being a 
measure for frequency instability, the ADEV is a well 
established tool to derive information about the type of 
noise that is present in the measurement, e.g. [3, 4, 5, 
6]. The ADEV shows a averaging behaviour for 

white frequency noise, i.e. for uncorrelated frequency 
values, and a steeper power law for white phase 
noise, due to correlations between the according 
frequency values in the latter case. 

The question if an estimate of the ADEV (Jylr = 
T) at the total length T of a measurement can be used 
for the determination of the uncertainty u has been 
addressed in [7]. There it is shown that the ADEV is 
equal to u for white frequency noise, while for white 
phase noise a constant correction factor of order unity 
has to be applied. It is thus justified to profit from the 
fast reduction of (Jyij) with in the case of white 
phase noise. 

However, the H-counter approach is not the 
only choice to measure frequencies, and in fact not 
always the optimum one. Especially the fact that H- 
counters only use two phase measurements - essentially 
detecting the times of the signal’s zero-crossings at 
the beginning and at the end of the measurement 
interval - makes this technique very susceptible to wide 
bandwidth phase noise. 

Eor this reason, counters have been developed that 
average the phase over parts of the interval, in order to 
effectively reduce the bandwidth and thus the influence 
of phase noise. In the extreme case, for an averaging 
interval with length r (actually data over an interval of 


at least 2t is required), the phase is averaged between 
times 0 ... T and r ... 2r, which would be called a A- 
counter. As we will show, the appropriate method 
to determine uncertainties associated to A-weighted 
averages is the modified Allan deviation mod Gy (r). 

Practically, modCTj,(r) has been introduced to 
resolve the ambiguity that is present in the functional 
dependence of the ADEV on r between white and 
flicker phase noise and to avoid a dependence of 
the ADEV on the measurement bandwidth in these 
cases [8]. As the modADEV shows an even steeper 
T“^/^-power law for white phase noise than the ADEV, 
the question arises if A-averaging could be employed 
advantageously to achieve a reduced uncertainty u by 
using the associated modADEV for its determination. 

In many cases, white frequency noise dominates 
at long averaging times, and thus a rejection of phase 
noise by the modADEV is not necessary [5]. As 
we will show, the use of the modADEV is even 
slightly detrimental in this case. There are, however, 
important measurement examples, for which white 
phase noise dominates even for r = 10^ s... 10® s [9, 
10, 11, 12], and the uncertainty u can be reduced 
significantly by a better data averaging and uncertainty 
analysis strategy, as for example by A-averaging in 
conjunction with uncertainty determination employing 
the modADEV. 

In section 2, the relationships between the 
uncertainty of H- and A-weighted averages of a 
measurement data set and the ADEV and modADEV 
are derived, since these are the technically widely 
implemented types. In section 3, the conversion 
factors from the modADEV to the uncertainty of 
the corresponding A-weighted average are deduced for 
different noise types. The question of an optimum 
averaging method is addressed in section 4, whereas the 
topic is discussed in terms of correlations in section 5. 
In section 6, a strategy for the transition between white 
phase and white frequency noise is povided. Finally, 
some practical issues are discussed in section 7. 

2. Average frequency and associated 
uncertainty 

The value of an oscillatory signal with amplitude Vq) 
phase d>(t) and nominal frequency t'o is described by 

V(t) = Vo sin [4>(t)] = Vo sin [2TTVot + (j){t)]. (1) 

The residual phase (j){t) contains the fluctuations due 
to noise and thus the statistical properties of the 
signal, whereas the nominal frequency i>o is defined as 
a time-independent value near the “true average” of 
the oscillator frequency. This true frequency cannot 
be perfectly known as a matter of principle, but 
only be estimated by performing a finite number of 
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experiments. The instantaneous frequency is defined 
by 


v{f) = 


1 d$ 
2ti dt 




1 d(j) 
2 ti dt 


( 2 ) 


Usually, the fractional frequency deviation y{t) is 
defined as 


Vo 


1 A(j) 

2ttvq dt 


( 3 ) 


The targeted result of a frequency measurement is an 
estimator y for the expectation value /i = (y) of the 
fractional frequency, along with a suitable uncertainty 
u{y). The expectation value is defined as the average 
over all elements of a fictitious ensemble of an infinite 
number of independent realizations of the random 
variable y. Throughout this paper, the angled brackets 
indicate the ensemble average. The uncertainty is 
defined as the scatter of the weighted averages y^ 
around the expectation value y, where the y^ are the 
realizations mentioned above. This scatter is described 
by the classical standard deviation, i.e. the squared 
uncertainty is defined by the classical variance 

= {iVw - h-f) = {iVw - {yw)f) ■ ( 4 ) 


As mentioned in the introduction, this “statistical” 
uncertainty does not include other contributions to 
the overall uncertainty besides those from noise. To 
determine u{y), the classical variance can be estimated 
by analysing the noise properties of the underlying 
measurement data. In the frequency domain, the 
power spectral density Sy{f) contains the necessary 
information about the fluctuations, whereas in the 
time domain the Allan variance cFy{T) and its derived 
variants are usually employed to analyse the noise 
properties. 

For convenience of reading, we will from here on 
not distinguish explicitly between the true uncertainty 
given by the classical variance and its estimate. This 
does not change the presented results for practical 
purposes. Furthermore, we presuppose without 
further discussion at several steps of the mathematical 
derivations that the underlying noise processes are 
ergodic. 

In the following, we will first treat the estimation 
of the expectation value y using 11- and A-weighting 
as prominent examples of averaging methods. In a 
second step, we derive equations relating the power 
spectral density with the uncertainty u associated to 
the averaging method. In a third step, we derive 
equations relating the power spectral density with the 
Allan variance variants corresponding to the respective 
averaging methods. Finally, we compare the equations 
obtained in the second and third step, in order to find 
the noise conditions under which the uncertainty u 
can be expressed using the according Allan deviation. 


and to find the appropriate conversion factors between 
them. 

We describe an averaging method by means of 
a weighting function w{t), which is non-zero only 
over a finite time interval between t = —t/ 2 and 
t = T12. Furthermore, we define a continuous 
moving weighted average function (t) given by 
the convolution between the instantaneous fractional 
frequency y{t) and the weighting function. Vwit) 
describes a sample of y{t) averaged over an interval 
of length T around t 

/ OO 

y{t')w{t'- t) dt'. (5) 

-OO 


We see from equations 4 and 5 that u{y) is intrinsically 
connected to the weighting function, i.e. to the 
method used for the determination of the average 
value. Different weighting functions are known, e.g. 
for rectangularly weighted, i.e Il-averaging as with an 
ideal Il-counter 


wn{t) 


Ijr —T12 < t < T12 
0 elsewhere 


( 6 ) 


and for A-averaging 

r I -I- —r < t < 0 

Whit) = \ 1 — tfr'^ 0 < t < r 
I 0 elsewhere 


( 7 ) 


Another interesting averaging method is to perform a 
least-squares fit of a linear function to the phase /fit) 
over the interval —t/2...t/2 [13, 6]. According to 
equation 3, the average fractional frequency at t = 0 
is given by the slope of the regression line divided by 
27rr'o- For the regression, 

= [ [4>{t) - (</o + ‘2T^voy^t)f dt (8) 

J-t/2 

is minimized with respect to the parameters (fo and y^ 

o 2 

of the regression line. Hence, = 0 and 


dx 


I'T/Z 

-Attvq / tcfit) - tcfo - 2nvoy.^f dt = 0, (9) 
J~t/2 

which implies 


12 


.r/2 


2^121 


0 J-t/2 


2TTVn 


t(j){f) dt 


r/2 


dt 




f- t/2 


2 dt 


.( 10 ) 


where we have integrated by parts. Hence with 
equation 3 

”/2 / .q et^' 




-t/2 

.t/2 

/-t/2 


3 


y{f) dt 


Wps.T{t)y{t) dt 


( 11 ) 
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and the frequency weighting function for linear 
regression of phase data is found as 

3/(2r) — < t < t 12 

0 elsewhere 


WQ.(t) = 


( 12 ) 


In analogy to the well-known 11- and A-counters and 
corresponding weighting functions, the nomenclature 
fl-counter and f2-weighting function is used due to the 
parabolic shape of the weighting function associated 
with a linear regression to phase data [14]. The 
temporal weighting functions are shown in figure 1. 

Plugging equation 5 into equation 4 yields for the 
uncertainty 

(13) 



Using Parseval’s theorem to express this by the single¬ 
sided power spectral density and Fourier transform 


W{f ) = J^^exp{2nift)w{t) dt of the weighting 


function yields 


= / Sy{f) \W{f)\^ df 

(14) 

with 



(15) 


(16) 

and 


TiA 3sin(7r/T) 3 cos (nfr) 

(17) 


The Allan variance is the two-sample variance without 
dead time 

= \{{yw„{t + T/2)-y^^{t-T/2)f'^ , (18) 

where the two samples are the fractional frequency 
averages taken over two intervals of duration r 
separated by r, i.e. without dead time. We now develop 
the mathematical justification and noise conditions 
for determining the uncertainty from the two-sample 
variance associated to a frequency averaging method 
described by a weighting function w(t). 

By inserting equations 5 and 6 into equation 18, 
the ADEV can be expressed as 

(19) 


= 


y{t')w„{t' -t) dt' 


with the weighting function for the ADEV given by 

Wa{t) = ^ [wn(t -f t/2) - wu{t - r/2)]. (20) 

The definition of the ADEV in equation 18 shows the 
connection between the ADEV and Il-averaging: The 
Allan variance is defined as the two-sample variance of 



Figure 1. Weighting functions for If- (red), A- (green) and £7- 
(blue) counters (left), and for the corresponding Allan deviation 
types (right). 


Il-averaged samples. This can be used to formulate 
a generalized definition of different types of Allan 
variances as the two-sample variance, where the two 
samples are determined using the same weighting 
function as it is used for the average value, whose 
uncertainty shall be determined. This generalized 
definition of the two-sample variance for a certain 
averaging type with xxx used as a placeholder for the 
type description reads as 


^ C'y ("t) = ^ ( (^ + 'r/2) - (t - t/2)] ^ ^ 


/ oo 

y(i')w^xxx<T(t' - t) dt' 

-oo 


( 21 ) 


with the one-sample weighting function Wxxx(i) and the 
two-sample weighting function w^xxo-if) associated to 
the averaging type xxx, given by 


ccr(i) [dix 


,(t-h t/2) - 'u;xxx(i -'r/2)]. (22) 


For example, if the averaging type is A, the associated 
two-sample variance is the modADEV, and for phase 
regression averaging it is the parabolic Allan deviation, 
abbreviated here as par ADEV [14]. 

Applying Parseval’s theorem to equation 21, the 
generalized Allan variance can be expressed by the 
single-sided power spectral density Sy{f) and Fourier 
transform W{f) of the temporal weighting function 
w{t) [15]. In the general case, the difference between 
the two samples time-shifted by t as they are used 
in the two-sample variance leads to an interference 
term in Fourier space, i.e. to multiplication of the 
according spectral weighting function W{f) by the 
factor •\/2sin(7r/r) 


pOO 

xxxct^(t)= / 5y(/) 1 Wxxx<t(/)|^ d/ 

^0 

poo 

= 2/ 5y(/)|Wxxx(/)sin(T/T)l" d/. (23) 

Jo 

From this equation the well known power laws for 
the scaling of modADEV (and similarly ADEV) with 
averaging time r can be calculated for different noise 
types and their corresponding power spectral densities 
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To establish the connection between the two- 
sample variance (e.g. modADEV or ADEV) and the 
uncertainty u(ij), equations 23 and 14 are compared 
with each other, leading to the values detailed in the 
next section. 

3. Application to different noise types 

To validate that the rapid fall-off of the modADEV and 
parADEV for white phase noise can be exploited by 
using A-averaging in conjunction with the associated 
modADEV /parADEV as an advantageous measure of 
uncertainty, the integrals given above were calculated. 
Eor completeness, we also provide the already known 
results for the ADEV. 

In the case of the ADEV, a high frequency cut¬ 
off fh ( 27 rT)“^ has to be introduced to ensure 
convergence of the integrals [8]. In practice, high 
frequency noise is always cut off by the limited 
bandwidth of the components used or by deliberate 
low pass filtering. Eurthermore, the ratio between uu 
and ay is independent of fh such that it does not affect 
the conclusions here. 

We have listed the results for white phase, flicker 
phase, white frequency, flicker frequency, and random 
walk frequency noise in table 1. Under the assumption 
fh we have neglected a small sinusoidal 

contribution depending on fh in the corresponding 
entries. It can be seen that for the first three 
noise types the power law of the respective Allan 
variance type is identical to that of the corresponding 
uncertainty u. We draw two main conclusions: 

(i) In the presence of dominating white or flicker 
phase noise one can indeed profit from the 
rejection of high frequency phase noise for the 
determination of the average y using A- or D- 
weighting . 

(ii) The modADEV / parADEV can be employed 
to determine the uncertainty, provided that A- 
/ f2-averaging and the appropriate correction 
factor listed in table 1 are applied. This is 
not conceptually different to using the ADEV to 
determine the uncertainty, but leads to tighter 
values of u if white phase noise dominates. 

In practice, the question arises how to estimate the 
according Allan deviation for the determination of the 
uncertainty u{y) of the average y taken over the entire 
available measurement length T. In the following 
we assume that the values of the Allan deviations at 
sampling times t = T and below are known, e.g. from 
other measurements, as it was also presupposed in [7]. 
Since the A-weighting function is nonzero over twice 
the averaging interval, the average value y and its 
uncertainty ua can be determined over an averaging 
interval of maximum length T /2. Instead, for the other 


Ratio 

White 

PM 

Flicker 

PM 

«/(T/2) 

4 

161n(2) 


fhT 

3 ln(27r/hT) —ln(2)+37 

«n(T) 

3 

9 


fhT 

3 ln(27r/),T) —ln(2)-|-37 

uliT) 

3 

9 ~ n yi o 


4 

161n(2) ~ U.OlZ 


White 

EM 


< 1 


4 

3 


< 1 


6 

5 


_9_ 

10 


Table 2. Ratios of the squared uncertainties obtained for A-, If- 
and tl-averaging using the respective available averaging times 
that can be obtained from a given dataset of length T, under the 
assumption that fh (27rr)“^. 


weighting functions, which are restricted to exactly the 
averaging interval, the uncertainty can be determined 
for averaging time T. 

The ratios u\{T / 2) / u^{T), Uq{T)/u‘^{T) and 
Uq{T)/u\{T/2) for the uncertainties determined by 
the described procedure are listed for typical noise 
types in table 2. Although for white and Elicker 
phase noise, some of the ratios depend on the 11- 
averaging cutoff-frequency fh, T is under practical 
conditions large enough to satisfy the conditions that 
A- and f2-averaging are advantageous over Il-averaging 
in these cases. However, this is not true if white 
frequency noise is dominating, where the application 
of A- or H-type averaging is actually detrimental 
(ratio > I), even though the direct comparison between 
ADEV, modADEV and parADEV in table I might 
not suggest that. This shows that it is important for 
the determination of the uncertainty to be aware of 
the difference between the classical variance and the 
two-sample variance such as the ADEV, modADEV or 
parADEV. 

Comparing fl- with A-averaging reveals that H- 
is in principle advantageous over A-averaging for 
all three noise types. However, in practice there 
are not yet many O frequency counters. Eor that 
reason, H-averaging at present does not yet play an 
important role, although this may change in future 
with increasing digital signal processing capabilities. 
For the choice of the averaging method in post¬ 
processing-based average determination, H-averaging 
is definitely a considerable option. 

4. Optimum averaging methods 

The methodology described so far can be readily 
transferred to any other averaging method, e.g. 
describing counters with non-zero dead times [6, 7] by 
adaptation of the weighting function. We now derive 
the optimum weighting functions for a given noise 
type by first recalling the general requirements and 
constraints imposed on the weighting function. (I) It 
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White PM 

Flicker PM 

White FM 

Flicker FM 

Random Walk FM 

Syif) 

h2p 

hif 

ho 

h-if-^ 

h-2f-^ 

<7y(E) 

U ^ — 2 

4772 “2” 

ln(27r/hT)-|-7 __2 

27r2 

^hoT-^ 

21n(2)/z_i 

|7r^/l_2T 

WnM 

^/i2r-2 

= hlir) 

31n(2ir/hT)-ln(2)-|-37 ^ __2 

\hoT-^ 

00 

CX) 

« hlir) 

modcr^(T) 


241n(2)-91n(3). __2 

\hoT-^ 

iln(2)h_i 

^n^h-2T 

uKt) 

4^h2T-^ 

= |modcry(r) 

« 0.822mod ct^(t) 

\hoT-^ 

= |mod(T^(T) 

00 

00 

par cr2(T) 

2f^^2r-3 

121n(2)-3. 2 

2772 “1^ 

|hoT-i 

14-8 ln(2)^ 

§TT^h_2T 


2f^^2r-3 
= par 0-2 (r) 

ss 0.846 par (t^(t) 

|hoT-i 
= par 0-2 (r) 

00 

00 


Table 1. Summary of the uncertainties u associated to If-, A-, and Q- weighted averages, in comparison with the ADEV (o'y), 
modADEV (modo-y) and par ADEV (par ay) as a function of the averaging time r, and for different noise types Sy{f), under the 
assumption ^ 7 Ri 0.577 is the Euler-Mascheroni constant. PM stands for phase noise, FM for frequency noise. 


leads to a minimum uncertainty, (2) must be confined 
to a finite averaging interval r in the time domain, and 
(3) must be normalized to 1 for integration over this 
finite length interval. To minimize the uncertainty, i.e. 
the classical variance for a given power spectral density 
Sy{f), we have to minimize the integral in equation 14. 
For white phase noise Sy{f) = h 2 p, we find 

pOO 

U^=h2 f\W{f)fdf 

Jo 

h 

= |27ri/lF(/)|" df 

= 

27r2 

= ^ L(f,w,w)dt, (24) 


dw{t) 


dt 


dt 


which has a functional form that can be minimized 
using the Euler-Lagrange formalism. We use the 
constraint 

.r/2 

/ w{t) dt = 1 (25) 

J-r/2 


as Lagrange multiplier to define the function 
Lx = L{t, w,w) — X (w - 

V T 

= — X (w — - 


(26) 


and write down the corresponding Euler-Lagrange 
equation 


dLx 

dw 


d dLx 


= 0 —A — 2w = 0. 


dt dw ^ ’ 

Integration of equation 27 using the boundary 
conditions w(< ± t/2) = 0 yields the weighting 
function given by equation 12. Hence, H-averaging 
yields the minimum uncertainty when white phase 
noise dominates. Independently of and simultaneously 
to this work, it was realized by Vernotte et al. that 
the parADEV is the optimum choice to determine 
the frequency instability in presence of white phase 
noise [16]. 

For white frequency noise instead, we can write 
equation 24 with L = jwp, and accordingly, 
equation 27 reads 


2r(; — A = 0 ^ w{t) = l/r for 0 < t < r. (28) 

Hence, the optimum weighting function in case of white 
frequency noise is the H-weighting function. 


5. Correlated frequency data 

It is instructive to investigate the subject from 
a slightly different point of view, that reflects a 
usual experimental situation: A counter with given 
averaging time r and chosen weighting function wn/A 
measures N consecutive frequency averages y^yi^lr), 
that shall in turn be averaged in a post-processing 


r 
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step to obtain the average y determined over the 
overall measurement time. The situation here is 
explicitly that the individual averaging windows are 
contiguously connected in the case of Il-averaging or 
even overlap with each other for A-averaging. This has 
direct consequences for the correlations between the 
individual averages 

The uncertainties u{y^) associated to the individ¬ 
ual averages are known from the above analysis or via 
a direct estimate of the classical variance (equation 4), 
and they can be assumed to be all equal in case of er- 
godic noise. The usual approach for the estimation of 
the combined uncertainty u{y) of the average y of the 
dataset consisting of N independent measurements is 
to multiply u{y^) by I/VN. 

However, in the general case correlations between 
the individual averages y^ilr) have to be investigated 
and taken care of [17, 7]. We will also see that for some 
situations, e.g. in which frequency data obtained with 
A-counters are subsequently H-averaged, this approach 
leads to results more easily than an adaptation of w{t) 
followed by the necessary integrations. The following 
method will also be used to optimize the averaging 
strategy for combined noise situations (section 6). 

According to the “propagation of uncertainties” 
e.g. described in the GUM [18], the combined 
uncertainty u{f) of a function f(xi) is found by 
linearization of / around the expectation values p-i and 
depends on the uncertainties of the N input quantities 
xi as 


N 


i(f\ dfdf 

(/) = 


l^m—1 
N 


dxi dXr^ 




dxt 


^-1 N Cif Cif 




dxi dxr, 


(29) 


with the covariance between xi and Xm 


uixi^ X-fri) — ((^/ yi^{Xm ^771')') • (^^) 

For the special case oi f = y = 
being the H-average of the xi = y^ilr) having identical 
uncertainties u{y^) and a common expectation value 
y, we find [17, 7] 


u^iV) 


1 

N 


U^iVw) 


2 


N-1 


^ {N-l)CyJlT) 

1^1 


U^iVw) 

N 


N-1 




CvJlr) 

CyJO) 


(31) 


using the autocovariance function Cy^ (t) of y^i,(t) — y, 
i.e. of the fluctuations of the fractional frequency y{t) 




H-counter 

A-counter 

Noise Type 

Syif) 

1 = 1 

l > 2 

1 = 1 

l > 2 

White PM 

h2p 

1 

2 

0 

1 

2 

0 

White FM 

ho 

0 

0 

1 

4 

0 


Table 3. Calculation of the autocorrelation function 
( 0)5 using equation 33 for different noise types, 
weighting functions and lag multipliers /. 


Noise Type 

H-counted 

A-counted 

White PM 

“n 

Pa 

Ar2 

White FM 

“n 

N 

3iV-l 2 ^ 3 

2Af2 “A ~ 2 Af 


Table 4. Overall average 'U^(y) for a set of N fl- or A-counted 
averages under the assumption of white phase or white 
frequency noise. 


convoluted with the weighting function w{t). It is 
defined as 


= i{yw{t + t)- n] [y^{t) - y]). (32) 


The Wiener-Khinchin theorem can be used to express 
the autocovariance function Cy^ (It) as the Fourier 
transform of the relevant double-sided spectral density 
%„(/) of Vniit) - M 


/ OO 

-OO 

^OO 

= / cos(2tt flT)sy^(f)df 


cos(2t fir) |Wn/A(/)|" Sy{f) df, (33) 

where in the second line we have used that the 
double-sided spectral density is even, because y^jit) 
is real. In the third line, the double-sided spectral 
density Sy^(f) has been substituted by the single¬ 
sided spectral density given by Sy{f) weighted by 
the spectral weighting function Wu/Aif) and the 
integration interval has been changed accordingly. 

Table 3 lists the values of the autocorrelation 
function defined by Cy^{lT)/Cy^(0) and computed 
from equation 33 for white frequency and white phase 
noise, H- and A-averaging, and lag It. We find that 
the autocorrelation vanishes for Z > 2 for the listed 
noise types (which is not the case for flicker phase 
noise) and that a A-counter introduces correlations due 
to overlapping averaging functions of the individual 
consecutive averaging processes, as can be seen in 
figure 1. 
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Inserting the results in table 3 into equation 31, 
we can determine the dependence of u^ify) on N. 
The respective results are listed in table 4. The 
dependence of u^ifj) on N expected for white phase 
and white frequency noise is retrieved for Il-counted 
frequency averages y^ [7, 17] and also found for 
the case of A-counters. We have thus built up a 
coherent theoretical framework that can be employed 
to find further optimized data analysis strategies in the 
presence of combined noise types. 

6. Optimization for a transition between noise 
types 

We have shown in section 2 that A-type averaging 
is advantageous over Il-averaging in the presence of 
phase noise to derive measurement results with small 
uncertainty u and that the modADEV can be used to 
determine u. Often, more than one type of noise is 
relevant in measurements. The combination of white 
phase and white frequency noise is a common example. 
Depending on the averaging time, one or the other 
noise type will dominate. Thus, according to table 1, 
using the advantageous uncertainty measure demands 
different average weighting functions in the r intervals 
differing in the dominating noise types. 

We will now develop an optimization strategy 
of mixed A- and Il-averaging to account for such a 
situation. For illustration it is helpful to consider 
a measurement situation that leads to an ADEV or 
modADEV as shown in figure 2 (left). A comparison 
of optical clocks via a long hbre frequency link [11] or 
via a satellite link can lead to similar plots. 

Calculating and plotting u (figure 2 (right)) 
according to section 2 for 11- and A-averaging 
illustrates the findings of tables 1 and 2. Due to 
the linearity of equation 14 in Sy{f), the uncertainty 
uffjyf) is the square root of the sum of squares of 
the uncertainties stemming from the individual noise 
contributions. 

Switching from A- to Il-averaging at an averaging 
time t' allows to unite the advantageous behaviour 
for both averaging approaches for the respective noise 
type (hgure 2 solid red line). Mathematically, this 
averaging approach can be described by an adapted 
weighting function wau or more easily by the approach 
outlined in section 5. We use A-averaged with 
r = t' = 10 s and average the N = T/t' — 1 data 
points that can be determined from an overall 
measurement length T. Explicit calculation of u^{y) 
in equation 31 for combined noise using equations 14 
and 33 reveals that ulfy) again is the square root of 
the sum of squares of the uncertainties due to the 
individual noise contributions. 


7. Practical implications 

The methodology described in the previous sections 
can be applied to frequency measurements by direct 
clock comparisons as well as to clock comparisons via 
some kind of link, e.g. a fibre or satellite link. The 
link-based comparisons are essentially also frequency 
measurements, but the drift and fluctuations of the 
link during a measurement enter the results as well. 

In the following, we discuss some more practical 
aspects occurring with real-world equipment used for 
frequency measurements, which are modern high- 
resolution frequency counters as well as satellite timing 
and ranging equipment (SATRE) modems. These 
aspects are closely related to further post-processing 
analysis of raw data generated by this equipment. 


7.1. State-of-the-art frequency counters 


In the past decades, frequency counter technology 
evolved [19, 5, 6] from conventional counting of integer 
periods during a certain gate time, via reciprocal and 
interpolating reciprocal counting, which adapts the 
gate time to an integer multiple of periods, to modern 
continuous phase recording based frequency counters, 
which are also referred to as continuous time-stamping 
frequency counters [20, 21, 13]. 

Essentially, the latter counter type measures 
frequencies obtained by dead-time-free Il-weighted 
frequency averaging over progressive time intervals of 
length To referenced to a signal with frequency /ref, i.e. 



(34) 


where fc is a fixed factor. The Il-weighted frequency 
average /n is achieved by a determination of the 
fractional number of cycles Cn of the signal under test 
that fits into the n-th averaging time interval, divided 
by the interval length tq. C„ is the phase difference 
accumulated over the averaging time interval tq in 
units of 27r. 


(fftn + r) - (f(tn) Cn Cn , 

/n c) 7 /refi 

Jttto To k 


(35) 


which shows that the phase recorder actually measures 
the frequency ratio between the signal under test 
and the reference signal. Furthermore, due to the 
absence of dead times between the averaging intervals, 
the phase values <j){tn) = J2n Cnj'T.'K are exactly 

given by the integration over the frequency values. 
The fractional number of cycles is determined by the 
combination of two techniques measuring the integer 
and the non-integer part separately. The integer part 
is determined with a continuously running counter, 
whose value is read out on-the-fly every tq. This 
strobed on-the-fly readout leads to the absence of dead 
times between the consecutive averaging intervals. 
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Figure 2. (left) ADEV and modADEV for data with a power spectral density Sy{f) of combined white frequency and white 
phase noise. Parameters are chosen such that the phase noise dominates for a significant part of the relevant measurement interval 
(ho = 0.01 s, /i 2 = lOs^, fh = 50 Hz, TO = 1 ms). For a discussion of the tq parameter see section 7.1. 

(right) Calculation of rt as a function of available data for the same noise spectrum as left using pure H-averaging (equation 6, dashed 
green line), pure A-averaging (equation 7, dashed-dotted blue line), and A-averaging up to a given r' = 10 s followed by H-averaging 
of 10 s A-counted values (solid red line). Using the optimized averaging approach, the white frequency noise limit (dotted black 
line) is approached at much shorter averaging times than for H-averaging. 


The non-integer part is measured with analogue 
interpolation electronics, e.g. based on charging a 
capacitor at a constant current [21]. 

More complex frequency averaging weighting 
functions for longer averaging times t > tq can be 
synthesized by proper addition and normalization of 
short If- weighting functions at wn(i+iiTo) of length tq, 
e.g. in order to approximate a A weighting function [5]. 
Alternatively, this can be done by according processing 
of the time series of the phase values as, e.g. by a linear 
regression for fl-averaging. 

Due to the granularity on the tq grid of 
the synthesized weighting function, the two-sample 
deviation associated to the averaging method (e.g. 
the modADEV associated to synthesized A-weighted 
frequency averaging) coincides at tq with the Allan 
deviation (t(to), i.e. with the two-sample variance 
associated to Il-weighted frequency averaging (see [5] 
for the case of a synthesized A weighting function). 
This value cr(To) strongly depends on the cutoff- 
frequency fh for the H-averaging. Hence, in order to 
make Allan deviation plots generated from different 
experiments comparable, it is important to provide 
the value of fh if an Allan deviation is derived from 
experimental data measured with H-counters, and to 
provide both the values fh and tq if another two-sample 
deviation (e.g. the modADEV) variant is determined 
using a “synthesized” weighting function for averaging. 

We consider as an example the K-l-K GmbH 
EXE frequency counter [21], which is a synchronously 
sampling multichannel progressive phase recording 
frequency counter nowadays used in many frequency 
metrology labs. It implements phase samples on a 


To « 1 ms grid. This value should not be confused 
with the resolution of the phase interpolation, which 
corresponds to the resolution at which tq is determined 
and has a value of 12.2 ps for the EXE. The EXE 
allows a large cutoff frequency fh of several tens 
of Megahertz. Thus, fh is usually determined by 
deliberate filtering of the signal before applying it 
to the input of the EXE. The EXE can save the 
raw phase values 4>{tn) time-stamped on the tq = 
1 ms grid to a file on a computer. Alternatively, 
the internal microprocessor of the EXE can perform 
an on-the-fly H- or A-averaging for averaging times 
choosable between 1 ms and 30 s and save them to 
a file on a computer for further post-processing. The 
synchronously sampled, in principle unlimited number 
of parallel counter channels of the EXE allows to 
measure differential phases or frequencies between the 
input signals. Fluctuations of the reference signal 
applied to the EXE are thus eliminated by a very high 
common mode rejection. 

1.2. SATRE modems 

As an alternative method to determine a coherent 
time series of phase values (j)„ on a tq grid, the 
sampled output of a phase comparator can be used. 
This is usually done in SATRE modems for two-way 
satellite time and frequency transfer (TWSTFT) [22]. 
It is straightforward that the output of the phase 
comparator depends on its nonlinearity, which may 
especially be critical for an analogue phase comparator. 
Nevertheless, noise can be filtered out either due to 
limited bandwidth of the delay locked loop employed 
for code locking or with a narrow low-pass filter behind 
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the phase comparator, and in order to obtain the phase 
samples, the temporal phase can be averaged according 
to 

_ I l-tn+T 

Htn) = — (36) 

h) Jtn —T 

which yields A-weighted fractional frequency averages 
given by 

— U \ 4‘{tn+l) ~ 4>{tn) 

V\{tn) = -Z-■ (37) 

27rTo 

The values of the y{tn) time series can be read out to 
a computer and written to a file as a time-stamped 
list on a To = 1 s grid. Again, these data can be 
further processed in order to obtain frequency averages 
with synthesized weighting functions. For example to 
approximate Il-weighted frequency averages by short A 
averages, we have to determine the difference between 
the last and first phase value of the (/>(t„) time series 
and divide by t. For a A-averaged frequency over the 
interval t, we have to subtract the equally weighted 
average over the 4>(tn) in the interval 0...T and 
over those in the interval —r... 0 from each other 
and divide by r. Please note that for A-averaging, 
there is no granularity involved in this case, because 
the elementary averages at tq are of the same kind. 
For linear regression weighting, the average angular 
frequency is given by the slope of the linear regression 
line fitted to the time series (j){tn) within the interval 
T. Regression weighting is particularly advantageous 
for clock comparisons via links with SATRE-modems, 
because it is robust against gaps in the phase data, 
which are always present using real satellite links [23]. 

8. Conclusion 

This paper provides the necessary mathematical 
foundation for the determination of the statistical 
uncertainty of a frequency measurement using the 
modified Allan deviation and other variants of 
the Allan deviation. When white phase noise is 
the dominating noise type, the steeper power law 
dependence on averaging time of the modified Allan 
deviation and of the parabolic Allan deviation with 
respect to the Allan deviation can indeed be exploited 
to achieve a tighter uncertainty at identical averaging 
times. This has direct benefits like, e.g., shorter 
measurement times required for comparison of optical 
clock frequencies using satellite links. For optical clock 
comparisons via fibre links [24, 10, 9], we see that A- 
averaging will allow to reduce the link instability to 
below the state-of-the-art instability of optical lattice 
clocks of 2 X 10“^®/-y/T/s after a few 10 s. For pure 
H-averaging however, several hours of averaging would 
be required. Keeping in mind that phase coherence 
has to be maintained over this interval, we find strongly 


relaxed requirements on the operational stability of the 
hbre links. 

Conversion factors between the modified Allan 
deviation and the uncertainty of the relative frequency 
have been found for several typical noise types. The 
presented formalism is not limited to the discussed 
noise types and the variants of the Allan deviation, 
but rather applicable to any two-sample variance 
associated to an averaging weighting function. 

We derived that for predominant white phase 
noise, linear regression weighting (fl frequency aver¬ 
aging) yields a minimum uncertainty, as does 11 aver¬ 
aging in case of white frequency noise. Furthermore, 
we have shown that employing the modified Allan de¬ 
viation for the determination of the uncertainty can 
also have benefits even with white phase noise on a 
white frequency background. A formalism to use an 
advantageous uncertainty measure in case of combined 
noise contributions has been presented, which includes 
the practically important situation of white phase noise 
combined with white frequency noise. 

Acknowledgments 

This work was performed within the framework of 
the Centre of Quantum Engineering and Space- 
Time Research (QUEST). Partial funding from the 
German Research Eoundation DEG within the geo-Q 
collaborative research center GRC 1128 and reasearch 
training group RTG 1729 is acknowledged, as well 
as funding within the ITOC and QESOCAS projects 
in the European Metrology Research Programme 
EMRP. The EMRP is jointly funded by the EMRP 
participating countries within EURAMET and the 
European Union. 

The authors are grateful for useful discussions 
with Burghard Lipphardt and Sebastian Raupach. 

References 

[1] D. W. Allan, Proc. IEEE 54, 221 (1966). 

[ 2 ] D. W. Allan, IEEE Trans. Instrum. Meas. IM-36, 646 

(1987). 

[3] J. A. Barnes et al, IEEE Trans. Instrum. Meas. IM-20, 

105 (1971), see corrections in NIST TN-1337 [25]. 

[4] P. Lesage and T. Ayi, IEEE Trans. Instrum. Meas. 33, 332 

(1984). 

[5] E. Rubiola, Rev. Sci. Instrum. 76, 054703 (2005). 

[o] S. T. Dawkins, J. J. McFerran, and A. N. Luiten, IEEE 

Trans. Ultrason. Ferroelectr. Ereq. Control 54, 918 

(2007). 

[7] W.-K. Lee, D.-H. Yu, C. Y. Park, and J. Mun, Metrologia 

47, 24 (2010). 

[8] D. W. Allan and J. Barnes, in Proceedings of the 35^^ 

Ann. Freq, Control Symposium (Electronic Industries 

Association, Ft. Monmouth, NJ 07703, 1981), pp. 470- 

475, for corrections see [25]. 

[9] S. M. F. Raupach, A. Koczwara, and G. Grosche, Brillouin 

amplification supports 1 X 10“^® accuracy in optical 



On the relation between uncertainties of weighted frequency averages and the various types of Allan deviations!! 


frequency transfer over 1400 km of underground fibre, 
arXiv:1504.01567 [physics.optics], 2015. 

[10] S. M. F. Raupach, A. Koczwara, and G. Grosche, Opt. 

Express 22, 26537 (2014). 

[11] K. Predehl et al.^ Science 336, 441 (2012). 

[12] A. Bauch et al.^ Metrologia 43, 109 (2006). 

[13] S. Johansson, in Proceedings of the 2005 IEEE Inter¬ 

national Frequency Control Symposium and Exposition 
(IEEE, New York, 2005), pp. 628-635. 

[14] E. Benkler et al, Agreement upon nomenclature for 

linear phase regression weighted averaging, 2015, At 
the IFCS/EFTF 2015 conference, the authors of this 
paper together with E. Rubiola and F. Vernotte agreed 
upon a unified nomenclature for the averaging by linear 
regression to phase data and the associated two-sample 
variance. The parabolic shape of the weighting function 
motivates using the terms Q-counter, fl-weighting 
function, and the terms parabolic Allan variance / 
deviation, par cr^ / par cr^, parAVAR / parADEV and 
PVAR / PDEV for the associated two-sample variance 
/ deviation. In [16], the terms quadratic variance and 
QVAR have been used before this agreement was made. 

[15] W. Riley, NIST special publication 1065, NIST (unpub¬ 

lished). 

[16] F. Vernotte, M. Lenczner, P.-Y. Bourgeois, and E. Rubiola, 

to be published in the Proceedings of IFCS/EFTF 
2015 special issue, IEEE Transactions on Ultrasonics, 
Ferroelectrics, and Frequency Control (IEEE, New York, 
2015). 

[17] T. J. Witt, Eur. Phys. J. Special Topics 172, 137 (2009). 

[18] Evaluation of measurement data - Guide to the expression 

of uncertainty in measurement, JCGM 100:2008 (2008). 

[19] J. Kalisz, Metrologia 41, 17 (2004). 

[20] G. Kramer and W. Klische, in Proceedings of the 2001 IEEE 

International Frequency Control Symposium and PDA 
Exhibition (SPIE, Bellingham, WA 98225 USA, 2001), 
pp. 144-151. 

[21] G. Kramer and W. Klische, in Proceedings of the 18th 

European Frequency and Time Forum, Guildford, UK 
(lET, London, UK, 2004), pp. 595-602. 

[22] W. Schafer, A. Pawlitzki, and T. Kuhn, in 31st Annual 

Precise Time and Time Interval (PTTI) Meeting (PTTI, 
Dana Point, GA 92629, 1999). 

[23] I. Sesia and P. Tavella, Metrologia 45, S134 (2008). 

[24] S. Droste et al, Phys. Rev. Lett. Ill, 110801 (2013). 

[25] D. Sullivan, D. Allan, D. Howe, and F. Walls, NIST 

tech, note 1337, NIST, U.S Department of Gommerce, 
National Institute of Standards and Technology (un¬ 
published), online available at http://tf.nist.gov/cgi- 
bin/showpubs.pl. 



